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We show how to decompose any density matrix of the sim- 
plest binary composite systems, whether separable or not, in 
terms of only product vectors. We determine for all cases the 
minimal number of product vectors needed for such a decom- 
position. Separable states correspond to mixing from one to 
four pure product states. Inseparable states can be described 
as pseudomixtures of four or five pure product states, and 
can be made separable by mixing them with one or two pure 
product states. 
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Entanglement, inseparability and nonlocality are some 
of the most genuine quantum concepts. While for pure 
states it is well established since long ago that the non- 
local character of the composite system is revealed in 
different but equivalent ways, the situation is drastically 
different for mixed states. For example, for pure states 
the violation of some kind of Bell inequalities H , or the 
demonstration that no local hidden variable models can 
account for the correlations between the observables in 
each subsystem, are equivalent definitions of non-locality 
||. But for mixed states, described by density matri- 
ces, such equivalences fade away. Consider a composite 
quantum system described by a density matrix p in the 
Hilbert space H a ®Ttb- In the frame set by the concepts of 
our starting sentence, product or factorizable states are 
the simplest possible. They are of the form p p = p a ® p&, 
i.e. for them, and only for them, the description of the 
two isolated subsystems is equivalent to the description 
of the composite system. Recalling that subsystems are 
described by the reduced density matrices obtained via 
partial tracing: p a — Tr b p ( p b = Tr a p), a density matrix 
corresponds to a product or factorizable state iff 



p = Tr b p <g) Tr a p 



P = Pp- 



(1) 



Also their index of correlation (or mutual information) 
defined in terms of von Neumann entropies of the system 
and subsystems, 



I c = Trp In p - Trp a In p a - Trp b In p b , 



(2) 



vanishes, and this happens only for them |3j]. Their 
subsystems are uncorrelated. Any state which is not a 
product state presents some kind of correlation. They 
are called correlated states. Quantum mechanics has 
taught us that there is a hierarchy of correlations, and 
the physics in the different ranks is different. The sim- 
plest correlated states are the classically correlated ones. 



Separable states are either uncorrelated or classically cor- 
related. Their density matrices can always be written in 
the form: 



Ps = ^PiPa 



® pbf, 1 > Pi > 0; ^2 



P, 



= 1, 



(3) 



i.e. as a mixture of product states. Their characteriza- 
tion is notoriously difficult. Thus, given a density matrix 
which is known to describe a separable state only very re- 
cently algorithms for decomposing it according to Eq. (||) 
have been found ^J||; besides, the decomposition is not 
unique. In fact, only recently Peres and the Horodecki 
family HQ] have obtained a mathematical characteriza- 
tion of these states, at least when the dimension of the 
composite Hilbert space is 2 x 2 or 2 x 3. For these cases 
the necessary and sufficient condition for separability is 
that the matrix obtained by partially transposing the 
density matrix p is still a density matrix, i.e. with only 
non-negative eigenvalues 



(p T «y >o 



P = Ps 



(4) 



For composite systems described by Hilbert spaces of 
higher dimensions, the positivity condition of p Tb is only 
a necessary one for separability |t]] . Following the hierar- 
chy of correlations, we find the states that are no longer 
separable, i.e. p ^ p s . These states are called "EPR- 
states" ||, "inseparable", "non-local", and sometimes 
"entangled" or simply "quantum-correlated" to empha- 
size that their correlations are not strictly classical any- 
more, though often these labels do not refer to exactly 
the same states. This confusion reflects the need of a fur- 
ther subclassification of the inseparable states according 
to whether they admit local hidden variables, whether 
they violate some kind of Bell inequality etc.. 

The issue we want to address here is whether any state, 
even if non-local, allows for some kind of local descrip- 
tion. We will see that this leads to novel physical perspec- 
tives about non-locality. Thus the aim of this contribu- 
tion is to decompose any separable or inseparable density 
matrix of a binary composite system of dimension 2x2 
in terms of only product vectors, and to give for all cases 
the minimal number of product vectors needed. In other 
words, we give the minimal local description of any state, 
be it separable or not. (Here and in what follows "local" 
means that it refers to the subsystems). More specif- 
ically, we will start proving that any separable density 
matrix can always be written as: 



1 



Ps = ^p l {\e i ){e l \ <8> \fi){fi\); 1 > Pi 



> 0: 



Yl Pi 

1=1 



= 1, 



(5) 



with 1 < n < 4, and we will determine the minimal n as a 
function of p s . This introductory result completes the re- 
sult n < 5 of Q and reproduces the result n < 4 of ||] in a 
completely independent way. Calling statistical mixtures 
of pure product states |e^} <8> e^) G H a , \fi) eWk, 
local mixtures, and calling the smallest n its cardinality, 
Eq. (^) says that any separable density matrix is a local 
mixture of cardinality smaller than five. We then come 
to our main results. First, any pure inseparable state 
(Pq = Pq) can be written as: 

Pq = (l + qi + q 2 )pi +) 

2 

- y^gi(|gi)(gi| ® l^iX^il); < ft < 00, (6) 
i=l 

with separable of cardinality 3. The subscript q 
means inseparable or quantum correlated. And, second, 
any non-pure inseparable state (p q > p 2 ) can be written 



Pq 



= (1 + q)pM - q(\g)(g\ <g> \h)(h\); < q < oo, (7) 



with ps separable of cardinality 3 or 4. We finally de- 
termine the cardinality of ps as a function of p q . As a 
consequence of our results any inseparable density matrix 
can be written as what we call a pseudomixture: 



p g =(l + q)p { s +) -qp\ >;0<q<oo 



(-). 



(8) 



of cardinality n = rt( + ) + n^ + ' and > being the 

cardinalities of p s + ^ and pi \ Then, in a nutshell, our 
main result is to determine for any state its representa- 
tion in form of a local (pseudo)mixture of minimal 
and then minimal n^ + \ Local pseudomixtures have an 
interesting physical interpretation. Eq.(^) for instance, 
says that any inseparable mixed state can be made sepa- 
rable by mixing it with some pure product state or, that 
its quantum correlations can be completely washed out 
with only one single local mixing preparation. 

Before proving all this, let us mention that local pseu- 
domixtures lead immediately to an unambiguous mea- 
sure of entanglement, 



E(p q ) = min q, 



(9) 



where ois defined in Eq.(g). It is unambiguous because 
in Eq. (]8[) only product states appear and thus E(p q ) just 
represents the minimal local mixing needed to wash out 
all entanglement. Minimizing q is however different from 
minimizing and then n( + \ which is what we do here, 
and we postpone its study and comparison with other 
entanglement measures [p|,pLT|— pL3[ for the time being. 



In order to prove Eq.(||) we need the following theo- 
rems: 

Theoreml. For any plane V\ in C 2 ®C 2 defined by two 
product vectors \v\) and \v^), either all the states in this 
plane are product vectors, or there is no other product 
vector in it. 

Proof: With the help of SU(2) <E> SU(2) transformations, 
\v\) and 1^2) can always be expressed so that: 



Vi = a x 



■01 



cos A 
sin A 



cos B 
sin B 



(10) 



with < A,B < 7r/2; A and B not simultaneously van- 
ishing and a±,/3i € C. All vectors in V\ are product 
vector iff sin A sin B = 0. If sin A sin B ^ 0, then the 
only product vectors contained in Vi are the generators 
of the plane \v\) and \v2) ■ 

Corollary. If p has rank 2 and is separable it can always 
be expressed as a statistical mixture of two pure product 
states and thus p Tb is also of rank 2. 
It suffices to see that for any separable p of rank 2 its 
range, 7Z(p), is a plane of type Vi- If it only contains two 
product vectors, then necessarily p = p\v\){vi \ + (1 — 
p) 1^2} (t>2 1 for some < p < 1. In the case that all vectors 
in IZ(p) are product vectors then its spectral decompo- 
sition gives us immediately the desired decomposition. 
Since in any case: 



P = p|ei/i)(ei/i| + (1 -_p)|e 2 /2)(e 2 / 2 |, 



(11) 



it immediately follows that p Tb is also of rank 2. 
Theorem 2. Any plane V2 in C 2 ® C 2 contains at least 
one product vector. Some planes contain only one. 
Proof: Consider the plane V2 generated by two orthog- 
onal vectors. Again, with the help of SU(2) (g> SU(2) 
transformations, it can be expressed as 



V2 — a 2 




(12) 



with A,B,C e TZ and j,6,a2,02 S C. Assume that 
none of the generating vectors is a product vector, that 
is AB ^ and C 2 AB + jS ^ 0. Then a vector in T 2 is 
a product vector iff 

a\AB + a 2 p2C(B 2 - A 2 ) - (3 2 (C 2 AB + jS) = 0. (13) 

With the above restrictions on A, B,C, 7 and 8, there 
is always at least one non- vanishing solution (i.e. 0:2, /3 2 
such that a 2 ^ 2 ^ 0) of Eq.@, There is sometimes only 
one non- vanishing solution (see also PH). 

We can now outline our procedure for finding the de- 
composition of a separable state into four pure product 
states. We will first prove that five pure product states al- 
ways do, and then present the slightly more cumbersome 
proof of going from five to four pure product states. The 
algorithm consists in subtracting a projector onto a prod- 
uct vector from p s or pj b in such a way that r(p s )+r(pj b ) 



2 



diminishes at least in one unity (here r(p) means the rank 
of p). We then repeat the procedure till the desired de- 
composition is obtained. Consider the most general case, 
a separable state p s such that both itself and its partially 
transposed matrix are of rank 4: r(p s ) — r(pj b ) — 4. As 
we shall see all the other cases are subcases of this one. 
Define now 



Pip) = 



1 



■(Pa-p|ei,/i>(ei,/i|); 0<p<l, (14) 



and 



l-p 



ipj b -P|ei,/r)(ei,/r|); < p < 1, (15) 



where |ei) £ Ti. a and £ Ti-b are completely arbitrary 
states. For p small enough both p and p Tb , are positive, 
and therefore, due to Eq.(|4j), separable. Let us denote 
by pi the smallest value for which a zero eigenvalue ap- 
pears in p(p) or p{p) Tb . Let us assume that for pi one 
eigenvalue of p(p) is equal to zero, i.e. r{p{p\)) = 3 and 
r(p(pi) Th ) = 4 (the same argument holds for the opposite 
case). Consider now a new product vector belonging to 
the range of p{pi), |e2,/2) € Tt{p{pi)), and define a new 
density matrix: 



pip) 



1 



l-p 



(p(pi)-p|e 2 ,/2)(e 2 ,/2|); 0<p<l. (16) 



As before, for p small enough, both p(p) and p(p) Tb are 
non-negative and thus separable. Let us denote by P2 the 
smallest value of p for which either p(p) or p(p) Tb develop 
a new vanishing eigenvalue. It cannot be p(p) unless, be- 
cause of the corollary, p{p) Tb develops simultaneously two 
vanishing eigenvalues. Therefore, it is in general p~ip) Tb 
which will develop a new vanishing eigenvalue, so that 



riPip2))=r{p{p 2 ) Tb ) = 'i. 



(17) 



As p{p2) has a decomposition of the type of Eq.(||) with at 
least three terms, and pip2) Tb has the corresponding par- 
tially transposed one, there always exists a product state 
satisfying: |e 3 ,/ 3 ) E 1Z(p(p 2 )) and |e 3 ,/ 3 *) £ K(p(p 2 ) Tb ) 
p>|,|lG|. Define now: 



pi?) 



1 



l-p 



(pip2)-p\e 3 ,f 3 )(e 3 J 3 \); < p < 1. (18) 



It is clear from the corollary that a p 3 exists such that: 

r(pip 3 ) > 0) = r(pi P3 ) n > 0) = 2, (19) 

and then it immediately follows that: 

PiPn) =P4|e 4 ,/4)(e4,/4 

+ (l-P4)|e B ,/5>(e B ,/B|; 0<p 4 <l, (20) 

completing thus the decomposition of any separable 
state. Therefore : 



p s = piPi + p 2 (l - Pl)P2 
+ Ps(l ~ P2)(l - Pi)P 3 

+ P 4 (l-P3)(l-P2)(l-Pl)f4 

+ (l-p 4 )(l-p 3 ){l-p 2 )(l- Pl )P 5 , 



(21) 



where Pi = |e^ , fi)(ei, fi\ are projectors onto pure prod- 
uct vectors. This proves Eq. (|) with n < 5. Notice that 
if r(p s ) + r(p^ b ) < 8 then n < 5. 

Let us now show that even when r(p s ) + r(p^ b ) = 8 one 
can always find a decomposition into four pure product 
states instead of five. To do this, we shall prove that 
there exists always at least one projector P = |e, /)(e, / 1 
and its partially transposed P Tb = |e, f*)(e, f*\ that can 
be subtracted from p s and p^ b respectively in such a way 
that positivity is preserved and the rank of both matrices 
diminishes simultaneously by one unit. Let us proceed 
defining as in Eq. (|l4| ), but for each of the five product 
projectors of Eq. (|2l|), the following five matrices: 



Piip) = iPs -P\ei,fi){ei,fi\); < p < 1; i = 1, ...,5. 

l-p 

(22) 

We will fix two sets of five values of p by the ten condi- 
tions 



r(pi(p = Si) > 0) = 3 
ripj b ip = s l )>0) = 3. 



(23) 



These conditions determine the maximal weights con- 
sistent with positivity with which the projectors Pi — 
fi)(ei, fi\ and Pj b = |e,, /*)(e l; f* \ can be subtracted 
from p 8 and p^ b respectively. We now show that it is 
impossible that Si < Si Vi or that Si > Si Vi. From (ilf l 
one knows the expressions for Si and si as defined above, 



1 



{ei,fi\p s x \ei,fi) 
1 



(24) 



If we call pi the probabilities for which Pi appears in p s 
(cf. Eq. (|)), then if, say, Si < Si Vi it immediately 
follows 



i=l i=l 

which from Eq. (|2j) reads: 

5 5 

^,Pi(ei> fi\P7 X \ei> fi) > J^P^-K^y 1 ^^' 
i=i i=i 

(26) 



or equivalently: 

Tr(psP s - 1 )>Tr(pJ b (pJ b )- 1 ), 



(27) 



3 



which cannot be. Thus at least for one i, say j, Sj > Sj . If 
they are equal, then subtracting this \ej, fj){&j, fj\ from 
p s in Eq. (H4) allows to reach 



r(p( Sj )>0)=r(p( Sj f b >0) = 3 



(28) 



in one step. If Sj > Sj then by connectivity of the space 
of product vectors and continuity of s and s as defined 



by Eq. (24) as functions of the states of this space, there 
exists one|e, f){e, f\ which has s = s and for which Eq. 
( p8| ) holds. Thus always a decomposition with four terms 
exists, and Eq. ([5]) has been proven with 



n = max(r(p s ), r(pj b )) < 4. 



(29) 



Let us now obtain our main results, which refer to 
inseparable states. From Eq. (Eh we know that 



infa(p Tb ) < 



Pq, 



(30) 



where c(p) means the spectrum of p. Let us prove that 
p Tb has only one negative eigenvalue. If there were two 
one could always find, according to theorem 2, a product 
vector |e,/} in the plane defined by the corresponding 
two eigenvectors, and for which obviously 



(e,f\pf\e,f) <0. 
But the above expression is equivalent to 
(e,/*|p ? |e,/*}<0, 



(31) 



(32) 



which is impossible, since p q > 0. We will call the eigen- 
vector of negative eigenvalue \N), i.e., 



p Tb \N) = -N\N); N>0. 



(33) 



We will now see that p q can be made separable by mixing 
it statistically with an adequate separable density matrix, 



pi \ i.e. 



1 



i + g 

where < q < 00 is such that 
1 



(Pq + Qp[ 



p{q) Tb = — q {pT b +qp { - )Tb )>0. 



(34) 



(35) 



We want to do this in a doubly minimal way. We want to 
choose pi to have minimal rank, and we then want to 
choose the minimal q, i.e. such that p{q) Tb just develops a 
vanishing eigenvalue {r(p(q) Tb ) < 4). Notice that due to 
the Hellmann-Feynman theorem |17[] the only eigenvalue 
of p^ b which can become zero by adding a non-negative 
operator is its negative eigenvalue. 

We will show how this is done as a function of the rank 
of p q : 

1. Assume r(p q ) = 1. Here p q represents an entan- 
gled pure state, which can always be written with the 



help of the SU(2) ® SU(2) transformations in its canon- 
ical form (cf. Eq. (|l|)) (e| = (cos A, 0,0, sin A) with 
cosAsinA > 0. It turns out that (N\ = ^=(0,1,-1,0), 

and that r(p q b ) — 4, as <r(p q b ) — {cos 2 A, sin A, 
cos A sin A, — cos A sin A(— —N)}. So, in this case, the 
minimal q satisfies r(p(q) Tb ) = 3. This implies that the 
rank of pi ^ cannot be one. Indeed, if it were one, as 
r(p q ) = 1 it would imply r(p(q)) = 2. But the two 
conditions r(p(q) Tb ) = 3 and r(p(q)) — 2 cannot be si- 
multaneously satisfied for a separable density matrix (cf. 
corollary). On the other hand a pi ' with r(pi ■*) = 2 
which does the job can always be found. It leads to 
r(p(q) Tb ) = r(p(q)) = 3. It can be implemented by choos- 
ing the two product vectors which statistically mixed rep- 
resent pi to be the vectors \gt, hi) given by the Schmidt 
decomposition of \N), \N) = Ci\gi,h*) + c 2 \g2,hl). This 
proves Eq. (||) with = p(q), q = q\ + 92 and where 
the result of Eq. ( p9| ) shows that the cardinality of is 
pi +] is 3. 

2. Assume r(p q ) — 2. Taking |e, /) £ TZ(Pq) which by 
theorem 2 always exists, we write p q in the form (l| 



-(|*><*|+p|e,/)(e,/|); p > 0, (36) 



Pq = 



where |^) is an entangled vector which belongs to lZ{p q ). 
Let us now prove that r(p^ b ) — 4. In order to do so 
write |^) in its canonical form |e). Consider the par- 
tially transposed of Eq. (pq). Recall (from the previous 
case) that |e)(e| b has three positive and one negative 
eigenvalues. The negative eigenvalue cannot be made 
to vanish adding the non-negative operator je, f*)(e, f*\ 
because then p^ b > 0, which from Eq. (|J) is incon- 
sistent with p q being inseparable. This, recalling that 
positive eigenvalues certainly cannot be made to vanish, 
proves r(p^ b ) — 4. This, in fact, always holds, so that 
r {plj b ) = 4 independently of r(p q ). It is now not too 
difficult to show that for any |e, /) always at least one 

pi ' = \g,h)(g,h\ exists which allows to satisfy Eq. ( |35| ) 
with r(p(q)) ~ r(p(q) Tb ) = 3. The upshot of this is that 

Eq. (Q) holds with ps = p{q) of cardinality 3. 

3. Assume r{p q ) = 3. As the previous case always 



allowed to find a pi ' with r(p ( s ') = 1 this is a fortiori 
true now too. This proves Eq. (^), but it is now not 

obvious whether it can always be done with a pi of 
cardinality 3. In fact, it cannot, as the analysis of the 
following counterexample shows: 



1 



Pq 



1+P1+P2 



(|*)(*| +Pi|ei,/i)(e 1 ,/ 1 | +P2|e2,/2)(e 2j / 2 |); p l : 

(37) 



with I*) - |e), (ei| - (fi\ - (1,0) and <e 2 | = (0,1). In- 
deed, none of the \g, h) vectors belonging to 7Z(p q ), which 
either have \g) — |e 2 ) or \h) — \ fi), does the job, and thus 



4 



r(p(q)) = 4. On the other hand it is easy to find exam- 
ples of p q for which r(p(qj) — r(p(q) Th ) — 3. Thus Eq. 

(Q) is proven but ps does not have always cardinality 
3. This parallels the ambiguity of n for separable states 
of rank 3, for which also sometimes n — 3 and sometimes 
n = 4. 

4. Finally assume r(p q ) = 4. In this case, obviously 
Eq. (0) holds for pg of cardinality 4. 

To summarize, we have proven that any separable state 
in C 2 ® C 2 is a local mixture of at most cardinality four, 
that any inseparable state in C 2 <g>C 2 is a local pseudomix- 
ture of cardinality four or five and that any inseparable 
state can be made separable by mixing it with only one 
single pure product state, except if it is pure, in which 
case it needs to be mixed with two pure product states. 
Therefore, when a state has only quantum correlations 
these can be made classical by mixing it with two pure 
product states, while when it has both classical and quan- 
tum correlations, mixing it with one single pure product 
state suffices to wash out all quantum correlations. 
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